Five new transit epochs of the exoplanet OGLE-TR-lllb 
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We report five new transit epochs of tlie extrasolar planet OGLE-TR-lllb, 
observed in tlie v-HIGH and Bessell I bands witli tlie FORSl and F0RS2 at 
tlie ESO Very Large Telescope, between April and May 2008. The new transits 
have been combined with all previously published transit data for this planet to 
provide a new Transit Timing Variations (TTVs) analysis of its orbit. We find 
no TTVs with amplitudes larger than 1.5 minute s over a 4- yea r obsery ation time 



baseline, in agreement with the recent result by lAdams et al.l fcOlOaf ). Dynam- 
ical simulations fully exclude the presence of additional planets in the system 
with masses greater than 1.3, 0.4 and 0.5 at the 3:2, 1:2, 2:1 resonances, re- 
spectively. We also place an upper limit of about 30 on the mass of potential 
second planets in the region between the 3:2 and 1:2 mean-motion resonances. 

Subject headings: exoplanets: general — transiting exoplanets: individual(OGLE- 
TR-111) 



Introduction 



The method of Transit Timing Variations (TTVs) has been proposed by several recent 



theoretical wor ks as of great potentia " 



etary svste ms (IMiralda-Escude et a" 



exomoons ( ISartoretti fc Schneider 



to d e tect additional exoplanets in transitin g exoplan- 
2002yAgoj_et al.ll2005l: [Holman et al1l2005h . and even 



19991 iKipping et al.l l2009a.bl ) . In a system with addi- 



tional planets, the predicted central time of a transit can change periodically by a significant 
amount of time if a second planet, the perturber, is in mean motion resonance with the tran- 
siting one. In the case of exomoons, TTVs of up to few minutes can be produced by moons 
larger than 1 depending on the physical parameters of the orbital system. Variations of 
other transit parameters, such as the depth or duration of the transit s, have been predicted 



to al so indicate the presence of additional planets in those systems (IMiralda-Escude et al. 
2OO2I ). Another application of the TTV technique is the detection of long-term orbital 



period secular variations produced by tidal interactions between the star and the planet, 
star oblateness and general relativity. Those effects are predicted to introduce changes in 
the orbital periods of the transiting planets with amplitudes of the order of 0.3-10 ms/yr 
dMiralda-Escude et al.l[2002[ Iflevl fc GladmanlboO?! : Ijordan fc Bakos al.l[2008l ) and would be 
detectable in 10-20 years timescales given the precision of the current techniques. 



^Based on observations made with ESO Telescopes at the Paranal Observatories under programme ID 
278.C-5022 
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These theoretical predictions have prompted the work of several observational groups 
who in the past few years have started to monitor various transiting exoplan ets from the 
groun d. More rece ntly, dedicated tran sit search space missions like CoRoT (IBarge et al. 



20081 ) and Kepler (IBorucki et all I2OIOI ) have also started to collect high duty cycle and 
long-ter m monitoring data which allow TTV studies of the new transiting planets they 



find (e.g. iBean II2009I ). Most of the findings from observational efforts are still preliminary, 



but include the potential detection of TTVs with amplitudes of up to 3 min for the hot 
Jupiter WASP-3b, which could be explained by the presence of a ~ 15 p lanet near the 
outer 2:1 mean motion resonance with WASP-3b flMaciejewski et al.ll2010l ). and the also 
preliminary detection of a —60 ± 15 ms/yr orbital period decay of the hot Jupiter OGLE- 
TR-113b by lAdams et al.l (j2010bl ). Other preliminary results include the hints of transit 
du ration and orbital incl ination variations in the hot Neptune Gliese 436b system reported 
by Coughlin et al.l (l2008h and the s hift of about 3 seconds in the orbital period of X0-2b 
reported by [Fernandez et al.l (120091 ). Seve ral other system , have been monitored without 
showing any clear evidence of TTVs, e.g. ISteffen fc Agol I (120051 ) found no eviden c e of a 
second planet using timing data of eleven transits of TrES-lb, and lAgol &: Steffen I (120071 ) 
ruled out Earth-mass planets in low order resonances in HD 209458b system. In a later 
work, iMiller-Ricci et al.ll2008al Jbl found no TTVs with amplitudes larger than 45 seconds on 
the transits of HD 209458b and HD 189733b based on MOST data, and lwinn et~aD J2009I ) 
found no variations from a constant orbital period for WASP-4b, after combining two precise 
mi d-transit tirn i ng me asurements with another fi ve ra easurements for this planet reported 
by iGillon et al.l ( 12009! ). Recently, iHolman et al.l ( 120101 ) announced the first unquestionable 
evidence of TTVs in the double transiting planetary system Kepler-9. 

OGLE- TR-lllb was the firs t hot Jupiter for which tentatively significant TTVs were 
reported by iMinniti et al.l (l2007l ). hereafter M07. The first tw o precision transit timing 
data points for this planet were published by IWinn et al.l (l2007l ). W07. Shortly after, M07 
pubhshed a third data point which deviated b y about 5 r ninute s from the predicted W07 
transit ephemerides. In a subsequent paper, iDfaz et al.l (120081 ) (D08) reported two new 
consecutive transits and combined their new data with the other three epochs to find TTVs 
with amplitudes of up to 2.5 minutes, which the authors suggested could be explained by 
the presence of a 1.0 perturbing planet in an outer e ccentric orbit to OG LE-TR-lllb. 
The most recent TTVs analysis publication on this planet (lAdams et al.ll2010al . AlO) reports 
six new transit observations and no TTVs for this planet with amplitudes larger than 71 ± 
67 seconds. 

In this new work we present five additional transits of OGLE-TR-lllb observed between 
April and May 2008, and we perform a new homogeneous timing analysis of all available 
16 epochs to further study the presence or absence of additional planets in this system. 
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In section [2] we present the new observations and the data reduction. Section [3] describes 
the modehng of the hght curves. In sections H] and |5] we describe the timing and possible 
parameters evolution. In E] we discuss the mass limits for a unseen perturber and finally in 
section [7] we present our conclusions. 



2. Observations and data reduction 

Between April and M ay 2008, we observed five transits of the exoplanet OGLE-TR-lllb 



with FORSl and F0RS2 flAppenzeller et al.l 119981 ) at the ESO Very Large Telescope. The 



first four transits were fully covered in phase. The fifth transit was only partially covered 
because the target reached the telescope's airmass limit, and is only complete between phases 
0.97 and 1.01, which includes the out-of transit baseline before the transit, the ingress and 
most of the bottom of the transit. The UT date of the mid-time of the transit, instrument, 
filter band, exposure time, airmass range and number of frames of each observation are 
summarized in Table [T] FORSl and F0RS2 are visual focal-reducer imagers composed of 
two 2048 X 4096 E2V/MIT CCD detectors mosaics with a pixel scale of 0.126 arcsec pixel"^ 
each (high resolution mode). The field of view (FoV) of each camera is therefore 4.25 x 4.25 
arcminutes, large enough to include OGLE-TR-111 and several comparison stars. FORSl 
and F0RS2 have the same wavelength coverage (3000-11000 A) but FORSl is optimized 
for blue wavelengths (< 5000 A) while F0RS2 is for the red (> 6000 A). Our first transit 
was observed with FORSl using a v-HIGH filter (Ae// = 557 nm) while the rest of the 
transits were observed with F0RS2 using a Bessell I filter (Ag// = 768 nm). A very close 
star appears partially blended with the target given the resolution of the instrument and 
the typical seeing conditions during the observations (Figured]). The center of the field was 
selected so that OGLE-TR-111 and several good comparisons stars would fall on a single 
detector, while also locating a bright nearby star out of the field. However, diffraction spikes 
from that bright star moved across the field of view, occasionally reaching the location of the 
target, as illustrated in Figure [TJ Our subsequent analysis revealed additional background 
noise in some of the images due to this effect. This problem was most evident in the 2008- 
05-03 light curve (Figure [2]) where pronounced bumps were visible in the bottom of the 
light curve of this night. Along the images obtained during this night, the diffraction spikes 
rotated trough the FoV reaching the comparison stars at different phases of the transit. The 
bumps also appeared in the light curves of these comparison stars. We use the peak of the 
larger bump in the light curve of the target to match in phase the bumps of all the other 
light curves (see Figure |2]). We average the light curves of the comparison stars in a small 
region where the bumps were more evident (between the horizontal lines in Figure [2]) and 
finally in order to remove the bumps, we substracted this average to the light curve of the 
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transit. In the other nights it was not possible to identify clearly the bumps produced by 
the diffraction spikes, and therefore we could not reproduce the process mentioned before 
but we attributed some of the noise in the lights curve to this effect. 

We worked with the processed data provided by the VLT pipeline which performs the 
bias and flatfield corrections. The times at the start of the exposure are recorded in the 
images headers, in particular we used the value of the Modified Julian Day keyword of each 
image and transformed it to Bary centric Julian Day (see section H] for details). 

Our photom etric analysis wa s done with the Difference Image Analysis Package (DI- 
APL) written by IWozniakI ( l2000l ) and recen tly modified by W . The package is an 



implementation of the method developed by lAlard fc Lupton and it is optimized to 



work with very crowded fields and/or blended stars, as is the case of OGLE-TR-111 (Fig- 
ured]). DIAPL models PSF variations along the X-Y coordinates, scales and subtracts the 
flux of the stars of a template image for each frame, among other calculations. One of the 
disadvantages of this code is the time required to complete the process, specially since it is 
not possible to verify the suitability of the input parameters and the quality of the products 
until the last steps (subtraction and/or photometry). In our case, due the large size of the 
frames and the large number of stars in the FoV, one iteration required up to 20 minutes. 
To eliminate a few nearby saturated stars in the field that hindered the photometry and to 
reduce considerably the processing time we worked on ~ 500 x 500 pixel subframes. Working 
with these subframes, DIAPL estimations of the PSF, background and flux levels are more 
representative of the vicinity of our target in each frame and in the reference image; noise 
as well as obvious systematics in the final photometry are considerably reduced. Reference 
frames for each night were constructed by combining the 20-27 best images (in terms of see- 
ing and signal-to-noise) depending on the night's conditions. We used aperture radii between 
4 to 10 pixels (in DIAPL's task phot, bash) to perform the relative photometry of the target 
and comparison stars (8 to 15 stars) on each subtracted frame. 

To obtain an abso l ute no rmalization, we performed aperture photometry (DAOPHOT 



/ ALLSTARS, IStetsonI (119871 )) in the reference frame of each night using a curve of growth 



analysis to select the ideal aperture radii. 

The resultant light curve contains some remaining systematic variations, which we have 
modeled using linear regression fits of the out-of-transit data points against the airmass, 
average FWHM of the point spread function and/or background level around our target of 
each frame. Doing this we were able to achieve an RMS precision of ~ 0.0013 — 0.0027 mag 
in the light curves (almost reaching the Poisson noise level in the best nights or doubling it 



^Thc package is available at 



littp://users.camk.edu.pl/pych/DIAPL/inclex.litinl 
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in the worst). 



3. Modeling Light Curves 

We fitted our five new light curves toget her with all the light curv es previously published 
by W07, DOS and AlO, using JKTEB OlB dSouth worth et al.l liool. The fit also includes 



the light curve by iPietrukowicz et al.l ( 120101 . PIO hereafter), obtained from the reanalysis of 
the VIMOS data published by MOT. PIO found a problem with the times reported by M07, 
which results in a mid-transit epoch time difference of ~ 5 minutes. We only consider the 
PIO analysis of this transit from this point onwards. 

Among the parameters fitted by JKTEBOP for each light curve are: the planet to star 
radii ratio (k), the inclination (i) and eccentricity (e) of the orbit, the out-of-transit baseline 
fiux (-Ffe), the mid-time of transit (Tc), the quadratic limb darkening coefficients {ui and ^2)? 
and the sum of the fractional radii, + r^. The terms and are defined as Tp — Rpl ^ 
and Ts = Rs/a, where Rp and Rg are the absolute stellar and planetary radii, and a is 
the orbital semi-ma jor axis . In th e cas e of the limb darkening coefficients, we fixed U2 to 
the values given by IClaret I ( I2OOOI ) and I Claret I ( 120041 ) for each observation's filter (for the 
v-HIGH filter we use the v coefficient) and only left Uir^ as free parameter during the fitting 
process described below, where x denotes the filter band. We performed the same fitting 
method using a linear limb-darkening law obtaining basically the same final Xred each 
transit, which reveals that the photometric precision of the light curves is not sufficient to 
distinguish between limb darkening laws. Also, to minimize potential degeneracies between 
parameters, we fixed the eccentricity and the longitude of the periastron of the orbit to zero, 
and the plan e t to s tar mass ratio to mp/rris = 0.00061, adopting the mass values derived by 
Santos et al.l tood ). 



JKTEBOP also allows for an statistical determination of the error of each parameter, 
as well as an analysis of the impact of systematics in the light curves, via Monte Carlo 
simulations. We ran 10^ Monte Carlo simulations adding random simulated gaussian noise 
to the input parameters to estimate the uncertainty of each parameter while also testing for 
potential correlations. When all the parameters described above are left free, there are clear 
correlations between k, i, and Uix, and also between Vp + rg and i, as illustrated in Figure |3l 
No significant correlation was observed between any of the other parameters. 



These correlations can be minimized by fitting the parameters in three steps, also run- 



http://www.astro.keele.ac.uk/ jkt/codcs/jktebop.html 
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ning 10^ Monte Carlo simulations on each step. First, we fixed the inclination of the orbit 
to the value obtained by AlO, i = 88.3° together with the corresponding value of ui^. from 
the Claret tables and fitted k and Vp + for each light curve. Next, we fixed k and + 
to the weighted average of the individual fits obtained in step I, and left i and Ui^ as a free 
parameters. Finally, we adopted the weighted average of the resulting inclinations and the 
corresponding uix for each filter (because we have only one transit observed with v-HIGH 
and one with V we adopted directly the results of JKTEBOP for its Mi^), and fitted only for 
T(.. In Tabled we summarized each step of the fitting process. An example of the histograms 
of the distribution of values for each parameter for the night 2008-05-12 are represented in 
Figure HI The adopted values of each parameter for the individual light curves are summa- 
rize in Table [31 The average values for the system based on all light curves are summarized 
in Table HI To test the consistency of the Monte Carlo er ror estimates, we compared it to 



the results of the prayer-bead method (iBouchy et al.ll2005l ). While the Monte Carlo method 



gives an idea of the white noise in the data, any level of red noise (time correlated noise, ) 
is best characterized by the prayer-bead method. The errors obtained by the prayer-bead 
method were, in general, larger than the Monte Carlo errors, showing that red noise is the 
dominant factor in the light curves. In Table [3] we show the ratio R of the errors estimated by 
the Prayed-Bead and the Monte Carlo method. We adopted as la errors of the parameters 
reported in Table [3], the larger values between these two error estimations. 

The Levenberg-Marquardt Monte Carlo (LMMC) fitting method implemented by JK- 
TEBOP, can present some disadvantages with respect to the Markov Chain Monte Carlo 
(MCMC), method used by several other recent TTVs studies. For example, LMMC can be 
trapped in a local r nimima or /and can underestimate the errors of the fitted parameters (e.g. 



Fisher et al.l l2009l : iDriscoU I l2006l ). Despite this, our results in Table [3] for the analysis of 
previously published light curves are fully consistent with the parameter fit values reported 
in Table 5 of AlO using MCMC. LMMC and MCMC are expected to yield similar results 
in well-behaved parameter space, i.e. with no multiple minima, as seem to be the case of 
our dataset. Even though, it is possible that the errors of LMMC best-fit parameters (pro- 
vided by the parameter distribution of the Monte Carlo iterations) can be under estimated, 
but our adopted parameter errors are dominated/scaled by the red noise contribution (see 
the R values in Table [3] ), indeed our errors are very conservative in comparison with AlO 
estimations. 
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4. Timing Analysis 

The mid-time for each transit derived in the previous section and listed in Table 2 was 
converted to Barycentric Julian Days, expressed in terrest rial time, i.e. Bjp(TT ), following 



the standard timing reference system recommendations by lEastman et al.l ( l2010l ). Since our 
analysis includes data from different instruments/telescopes, spanning over four years, and 
reduced by different groups we checke carefully for any possible systematics. As example, 
AlO has already pointed out how previously published results on OGLE-TR-lllb had not 
corrected the reported times for leap seconds (UTC to TAl) nor for the 31.184 second 
conversion between TAI and TT, where TAI is defined as the International Atomic Time. 
We have applied the same corrections to all the literature light curves used in our analysis. 

As an additional check to our reduction, light curve fitting and timing correction proce- 
dure, we compared our final BJD(TT) times to those published by AlO. Particularly valuable 
for this test is the 2008-05-12 transit epoch, which was independently observed by us and 
AlO. As illustrated in Figure O in spite of adopting completely different approaches to fit 
the light curves, all our derived transit mid-times agree well with the values obtained by AlO. 
When comparing the of our 2008-05-12 transit with AlO's, both times agree within the 
la error of our observation, although our mid-transit time occurs 124 seconds earlier. Four 
of the five transits we measure also produce mid-transit times on average 100 seconds earlier 
than the ephemerides predict (see Figure ED, although the differences are not statistically 
very significant. However, we have decided to further investigate the potential source of this 
discrepancy. First, we confirmed with the VLT staff that the times recorded in the image 
headers are the UTC times at the beginning of each exposure of our data (including all leap 
seconds). We further tested for any potential systematics between datasets by binning our 
data to 60 seconds, and removing the last ~ 50 points in our light curve to make it match 
the light curve sampling and phase coverage of AlO. The result was only a 10 second time 
shift in the resulting compared with the previous fit value. The two transits were ob- 
served in different filters, but we find no correlations between the limb darkening coefficients 
and Tc (see Figure [3]) that could account for this mid-time discrepancy. We attributed this 
difference to the red noise in the FORS light curve, possibly due to weather (the R value of 
Tc of this transit in Table [3] is almost twice the value of AlO light curve). 

Figure [7^ shows the updated Observed minus Calculated (0-C) diagram with all the 
final BJD(TT) mid-time values calculated from the literature light curves and our five new 
transits. The data show a linear trend which can be attributed to the accumulation of timing 
uncertainties with respect to the adopted transit ephemerides (D08). After removing that 
linear trend (Figure [7)d), the data are consistent with a constant period ephemeris equation 
of the form: 
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Tc = 2454092.80691 (25) [SJD] + 4.0144477(16) x N, (1) 

where T^. is the central time of a transit in epochs since the reference time Tq. This 
fit has a reduced of l-S? so the errors in Tc have been rescaled by a factor of \/x^ =1.18 
to make them consistent with = 1 (see errorbars in Figure [7)d), The new period is fully 
consistent with the one obtained by A 10. 



5. Analysis of additional parameters of the light curves 

We tested for possible variations of the physical parameters of the OGLE-TR-111 system 
using the values of the transit duration, T14, the inclination, the planet to star radius ratio, 
and the sum of the fractional radii derived for each transit with JKTEBOP. The value of 
each of those parameters over time is represented in Figure [HI 

There is no evidence of trends for any of the inspected parameters. However, it is 
noticeable that the results from the light curves observed by different groups appear clustered 
around the same values. We attribute that clustering to systematics introduced by the way 
the photometry is performed (i.e. aperture or PSF photometry, differential image analysis, 
and so on), and probably also by the way the light curve systematics are treated by the 
different groups. For example, D08 already pointed out that the depths of their transits 
were smaller than the average and that this was due to a reductio n artif act of the differential 



image subtraction techniques previously noticed by iGillon et all ( 120071 ). Although we used 
DIAPL instead of aperture photometry, our results are consistent with those of AlO. Notice, 
however, that systematics in the transit depths have no effect in the determination of the 
transit midtimes. 



6. Limits to additional planets 



Based on the timing constrain of our O — C diagram, i.e. no TTV variations with 
amplitudes larger that 1.5 minutes over a three year period, we run dynamical simulations 
to place limits on the mass and the semi- niajor axis of a po ssible orbital companion of the 
transiting planet using the mercury code (jChambers 1119991 ). The first step was to explore 
stable orbital regions by assuming a massless point particle over a range of initial semi-major 
axes, and fixing all the other input variables to the known physical parameters of the system. 
The orbital evolution of the massless particle was integrated over 10^ days. This test yields 
a strip of unstable orbits between 0.034 — 0.056 AU, where encounters between the test 
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particle and OGLE-TR-lllb would occur. For all the other orbits we calculate TTVs of 
the transiting body with the hypothetical coplanar perturber using a wide range of masses 
(0.1 < Mper < 5000 Mearth), Variable density (from Earth to Jupiter density depending 
on the mass) and semi-major axes (0.02 AU < a < 0.13 AU in steps of 0.005 AU) with 
~ 4500 simulations over 7 years. All the initial relative angles were fixed to zero. Near 
resonances, the steps in the variables were reduced to increase precision. Only the last 5 
years were used for the timing analysis, since that is about the same time span covered by the 
observations, and also to minimize any effects introduced by the choice of initial parameters. 
Then we calculated the central time of each transit during all the simulations. Similarly to 
what was done in section HI we did a linear fit of these central times and we defined the 
TTV of each simulation as the standard deviation of the central times with respect to this 
linear fit. We checked that the final period of the transiting planet did not change by more 
than 3a from its initial value due the gravitational interaction. With this method we were 
able to make a mass vs a diagram (Figure [9]-A) where the solid line represents TTVs of 
1.5 minutes. For comparison we also plot the mass limits of TTVs of 0.5 and 5 minutes 
(dotted and dash-point lines). T he dashed line corr esponds to the detectability limit placed 
by radial velocity observations (jSantos et al.l 120061 ). Our mass constrains are upper limits 
since for perturbers with an orbital eccentricity different from zero, the mass necessary to 
produce TTVs of the same amplitude will be lower. We confirm this by performing a set of 
simulations with e = 0.3 and using as input parameters the values {Mper, a) which produced 
TTVs ~ 1.5 min in the case e = (see Figure |9]-B). With this configuration the unstable 
region becomes wider due to encounters between the orbital bodies and the TTVs produced 
for a given mass were larger than in the case of e = 0. Same results were obtained setting 
the initial values of the longitude of the periastron different from zero (90°, 180° and 270°) 
and e = 0.3. Combining TTV RMS and radial velocities we can rule out the presence of a 
perturber body with mass greater than 1.3, 0.4 and 0.5 M^arth at the 3:2, 1:2, 2:1 resonances 
with OGLE-TR-lllb and lower the upper limit in the region exterior to the planet until 
~ 0.08 AU to companions of less than ~ 30 M®. 



7. Conclusions 

We present 5 new transit light curves of OGLE-TR-lllb. We homogeneously model 
all available light curves in the literature and search for any variation in the timing of the 
transits. With our updated ephemeris equation we find no TTVs with amplitudes larger 
than 1.5 minutes and therefore we rule out the presence of a companion in the 2:1, 3:2 and 
1:3 orbital resonances. If the system has an additional orbiting body, its mass has to be lower 
than 30 M^ if is located between 3:2 and 1:2 resonances. The mass limits we place with our 
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dynamical simulations based in the TTV data are lower than the limits obtained with radial 
velocities alone. We search for any trend in the duration, depth of the transit and inclination 
of the orbit but we do not see any clear evidence of variation with statistical significance. 
We point out that systematics of no evident source in the observations, reduction and/or 
analisys processes can induce differences in the values of the parameters obtained from the 
light curves and therefore a monitoring of transiting exoplanets carry out by the same group 
can contribute to reduce these differences. 
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Table 1. Observational information of each night. 



Transit Date 


Instrument 


Filter 


Integration Time [s] 


airmass range 


# images'^ 


2008-04-26 


FORSl 


v-HIGH 


30 


1.25 - 1.43 


323 


2008-04-30 


FORS2 


Bcssell I 


12 


1.25 - 1.59 


488 


2008-05-04 


FORS2 


Bessell I 


12 


1.25 - 1.75 


522(2) 


2008-05-12'' 


FORS2 


Bessell I 


4 


1.25 - 2.18 


601(94) 


2008-05-20'' 


FORS2 


Bessell I 


8 


1.29 - 2.07 


373 



''This transit was also observed by I Adams et al ] l|2010al 'l. 
''This transit has a incomplete phase coverage. 

'^The number of images doscartod in the analisys is shown in parenthesis. 
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Table 2. Summary of each step of the fitting process. See section [3] for details. 





Parameters 




Step 


Free Fixed 


Reference'' 


I 


/l , 7'p -j- T" s , jTc i 


AlO 




Ulx 


Claret 


II 


i, Ulx, Tc k, rp + ra 


Step I 


III 


Tc Tp -\- Ts 


Step I 




i, Ulx 


Step II 



''This column shows the origin of the adopted 
value of the fixed parameter. 



Table 3. Adjusted parameters for each transit using JKTEBOP code. 



Transit da.te'^ 


k 


R 




R 


i [°] 


R 


Mia; 


R 


Tc - 2450000 (BJD) 


R 


2 

^red 


2008-04-26(122) 


fl 1203('23'1 


0.6 


0926('17"l 


0.8 


88 35('14'l 


0.4 


81 ("05") 


0.5 


4582 56853('4S'l 


0.5 


0.17 


2008-04- rinfl2S"l 


1168fl5~) 


1.1 


0929('1 3") 


1.3 


88.22(12) 


1.3 


201' 08 "I 


1.3 


4586 5S2SSC72'1 


2.2 


0.57 


2008-05-04(124) 


0.1169(21) 


1.6 


0.0947(19) 


2.0 


88.49(15) 


1.5 


0.22(05) 


0.9 


4590.59853(94) 


3.0 


0.69 


2008-05-12(126) 


0.1151(29) 


1.7 


0.0927(19) 


1.6 


88.24(27) 


2.3 


0.00(28) 


3.0 


4598.6261(18) 


4.2 


0.66 


2008-05-20(128) 


0.1265(24) 


1.4 


0.109(17) 


1.2 


89.37(99) 


0.8 


0.13(13) 


2.0 


4606.65482(89) 


1.9 


1.34 


2005-04-09(-155) 


0.1179(87) 


1.1 


0.0990(69) 


1.2 


88.89(53) 


0.7 


0.89(08) 


1.5 


3470.5684(11) 


1.5 


1.08 


2006-02-21(-76) 


0.1266(15) 


0.8 


0.0934(08) 


0.8 


88.70(36) 


0.4 


0.45(05) 


0.9 


3787.70928(56) 


1.9 


1.15 


2006-03-05(-73) 


0.1246(12) 


0.8 


0.0934(08) 


1.2 


88.78(16) 


1.2 


0.61(04) 


1.1 


3799.75213(81) 


2.9 


1.08 


2006-12-19(-l) 


0.1122(19) 


0.8 


0.0899(12) 


1.0 


87.85(17) 


1.7 


0.03(25) 


2.4 


4088.7922(19) 


4.9 


1.3 


2006-12-23(0) 


0.1139(26) 


1.5 


0.0910(11) 


1.0 


87.97(21) 


2.1 


0.06(22) 


2.3 


4092.8056(15) 


3.7 


1.14 


2008-04-18(120) 


0.1245(31) 


2.1 


0.098(16) 


1.4 


87.92(81) 


17.3 


0.26(12) 


1.6 


4574.54272(83) 


2.3 


0.99 


2008-04-22(121) 


0.1190(21) 


0.8 


0.0944(16) 


1.0 


88.54(17) 


0.9 


0.40(07) 


0.8 


4578.55497(43) 


1.0 


0.96 


2008-05-12(126) 


0.1183(13) 


1.2 


0.0915(09) 


1.3 


88.14(10) 


1.4 


0.26(07) 


1.5 


4598.62754(47) 


1.9 


1.02 


2008-05-16(127) 


0.1214(16) 


1.5 


0.0922(14) 


1.8 


88.38(13) 


1.5 


0.35(07) 


1.7 


4602.64167(47) 


1.7 


0.99 


2009-02-17(196) 


0.1187(15) 


1.7 


0.0926(15) 


2.5 


88.27(17) 


3.0 


0.27(09) 


2.7 


4879.63863(55) 


2.7 


1.01 


2009-03-13(202) 


0.1236(09) 


1.0 


0.0893(06) 


0.9 


88.13(13) 


1.5 


0.28(08) 


1.4 


4903.72566(26) 


1.0 


1.01 



Note. — The values shown in parenthesis correspond to the errors in the last digits. 
^The Epoch number is shown in parenthesis. 

^R: Prayer-Bead and Monte Carlo errors ratio. See text for details. 
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Table 4. Final Values 



Parameter 


Adopted Value 


Error 


k 


0.1213 


±0.0004 


rp + Vg 


0.0917 


±0.0003 


i n 


88.29 


±0.04 


Period [days] 


4.0144477 


±0.0000019 


To{BJD) 


2454092.80691 


±0.00030 


U\v~HiaH 


0.81 


±0.05 


"1/ 


0.426 


±0.024 




0.313 


±0.033 




0.89 


±0.08 
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Fig. 1.— Portion of 0.5' x 0.5' images of the night 2008-05-03 observed with F0RS2 at 
VLT. The best image (FWHM~ 0.34") of the night is shown at the left panel and the worst 
image (FWHM ~ 0.55" ) is shown at the right panel. The location of OGLE-TR-111 is 
marked by the circle. Due to the seeing and pixel scale our target appears blended with a 
nearby star. Also one of the diffraction spikes of a very bright star is visible at the upper 
left corner of the best seeing image. Ocassionally, these spikes reached the location of the 
target contaminating the photometry. 
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Frame number 

Fig. 2. — Light curve of OGLE-TR-111 (top) and two comparison stars (middle) of the night 
2008-05-03 observed with F0RS2 at VLT. The contamination produced by the diffraction 
spikes of a very bright star is evident in the region enclosed by the solid vertical lines. The 
average flux of the comparison stars (bottom) after align the peaks of the bumps (vertical 
dashed line) was substracted from the light curve of the transit to remove the contamination. 
This procedure was only applied to this transit. 
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Fig. 3.— 3000 results of 10000 Monte Carlo iterations (using JTKEBOP with the data of 
the night 2008-05-12) which show the correlation between the parameters fitted from the 
light curves: k, Vp + Vg, Tc, i and the linear coefficient Uix of a quadratic limb darkening law, 
when all are left variable. The correlation coefficients between each variable are show in the 
bottom left corner of each box. 
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Fig. 4. — Histograms of the 10000 Monte Carlo iterations (using JTKEBOP with the data 
of the night 2008-05-12) for the fitted parameters on each of the : k, Vp + r^, Tc, i and the 
hnear coefficient uix of a quadratic hmb darkening law obtained after the steps I, II and III 
described in Section [3l The dashed lines show the fitted value and the ilo" errors which were 
compared with the red noise estimation using the prayer-bead method. The same analysis 
was performed on each of the 16 transits. The final value of k, Vp + rg, i, Uu' and Uu 
corresponds to the weighted average of these results, except for Ui^^high and Uiv where we 
adopted the results of JKTEBOP (see text). 
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Fig. 5. — Light curves of all transits of OGLE-TR-lllb. The solid lines show our best 
model fits produced by JKTEBOP. The filter, epoch number and author of the light curve 
is also indicated. 
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Fig. 6. — Comparison of the eleven mid-transit times, Tc, obtained by AlO and recomputed 
by us in this work. Sohd circles show the difference between the Tc measured by each 
group. The open circle shows the Tc for the new transit we observed in 2008-05-12 UT, and 
coincides with one of the transits measured by AlO. Although we find that the transit occurs 
124 seconds earlier, both results are consistent within the errors. 
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Fig. 7. — Panel-A: Observed minus Calculated diagram of the central ti mes of the transits of 



OGLE - TR-lllb. B l ack d o ts represent the tim es obtained with data fro m lPietrukowicz et al. 



(j2ninl ). IWinn et al.l fennTh . lPiaz et all tOO^ and lAdams et all (|2010aD, while the white dots 
are from the transits of this work. The dashed line represents a linear fit of the data. In the 
small box a zoom of the points between 110th and 120th epochs is shown. Panel-B: When 
the linear trend is removed no variations of more than 1.5 minutes are present. The errors 
were rescaled by a/x^ = 1.18 to make them consistent with a linear fit with xled ~ 
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Fig. 8. — Resulting values for the ratio of the radii (Panel a), sum of the fractional radii 
(Panel b), orbit's inclination (Panel c) and duration of the transit (Panel d), defined as 
the time between the first and f ourth contact, for each lig ht cu rve using JKTEBO P. Open 
and solid triangles correspond to iPietrukowicz et al.l (|2010[) andlWinn et al.l (120071 ) transits. 
Solid squares correspond to the tw o transits of Ipfaz et al. ( 2008 ) and solid dots represent 
the transits of I Adams et al.l (|2010aJ). Open diamonds correspond to the new transits of this 
work. 
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a-perturber (AU) 

Fig. 9. — Panel-A: Upper mass limits of an orbital perturber. These simulations were 
computed using e = 0. The solid line represents transit timing variations of 1.5 minutes. 
The dotted line and dash-point line represent TTVs of 0.5 and 5 minutes, respectively. The 
dashed line corresponds to the limits due the radial velocities observations. Vertical lines and 
gray strip indicate the orbital resonances locations and the instability region respectively. 
An orbital companion of OGLE-TR-lllb should have a mass in the region below the black 
solid line which corresponds to the mass limit imposed by the timing analysis. Panel-B: 
Transit Timing Variations with e = 0.3. If the eccentricity of the perturber represented by 
the solid line in Panel-A is increased, it will exhibit larger values than 1.5 minutes for its 
TTVs. Regions with TTVs below 1 minute correspond to unstable orbits with this new 
configuration. 



